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Recent developments in Quantum Chromodynamics (QCD) are reviewed on three major topics where nonper- 
turbative gluon excitations of the QCD vacuum and the physical properties of the confining flux play a central 
role: (1) New lattice results on the spectrum and wave functions of heavy QQgluon molecules, known as heavy 
hybrids, will be discussed. (2) Recent advances on the glueball spectrum in lattice QCD will be presented with 
some theoretical observations. (3) Progress in our understanding of the nonperturbative internal structure of the 
confining flux and its excitation spectrum will be reported. 



INTRODUCTION AND OUTLINE 

I will review here three topics with consider- 
able progress in phenomenologically relevant and 
theoretically challenging areas of lattice QCD. 
They are all related to nonperturbative gluon 
excitations of the QCD vacuum and the physi- 
cal properties of the confining flux. The limited 
scope of the talk does not include discussions on 
some exotic objects, like light hybrids ||,|^ and 
dibaryons or glueball phenomenology 
and some other topics which were reviewed, or 
presented elsewhere at the conference. The out- 
line of the talk is as follows: 

1. Heavy Hybrid Spectrum 

- Methods and Phenomenology 

- Born-Oppenheimer Approximation 

- NRQCD Hybrid Spectrum 

- Results from the MILC Collaboration 

2. Glueball Spectrum 

- Review of Recent Lattice Results 

- Composite Glueball Operators 

- Consituent Gluon Model 

3. The Confining Flux 

- Microscopic Vacuum Model 

- Casimir Scaling and Multiple Sources 

- Where Is the String Limit? 

- Flux Fission (String Breaking) 

The chromoelectric flux with static sources is 
the organizing theme of the selected topics. Its 
ground state energy explains the origin of the 
confining potential of heavy QQ pairs. The exci- 
tation spectrum of the confining (and confined) 



flux defines hybrid QQ potentials which bind 
quark- antiquark pairs into heavy hybrid QQgluon 
molecules. Higher excited flux states are un- 
stable against glueball emission requiring a de- 
tailed study of the glueball spectrum. The con- 
fined fiux, when stretched, is expected to reach 
the limiting behavior of an effective QCD string 
in geometric string variables which remains a 
topic of considerable theoretical interest. The 
stretched flux can also fission by quark-antiquark 
pair creation in the QCD vacuum with important 
phenomenological implications. Flux states with 
sources in higher color representations, or with 
multiple color sources, like three-quark and four- 
quark configurations, provide further information 
on fiux formation in the confining vacuum. 

1. HEAVY HYBRID SPECTRUM 

I will review recent results on the heavy hybrid 
spectrum from three complementary approaches, 
followed by a brief phenomenological note. 

Methods and Phenomenology 

The first method is based on the Born- 
Oppenheimer expansion. In this approach, the 
heavy hybrid meson is treated analogous to a di- 
atomic molecule: the slow heavy quarks corre- 
spond to the nuclei and the fast gluon field cor- 
responds to the electrons [^,1). The first step 
in the Born-Oppenheimer treatment is to de- 
termine the energy levels of the glue (and light 
quark-antiquark pairs) as a function of the heavy 
quark-antiquark separation, treating the heavy 



quark and antiquark simply as spatially-fixed 
color sources. Each such energy level defines an 
adiabatic potential. The quark motion is then re- 
stored by solving the nonrelativistic Schrodinger 
equation using these potentials. Conventional 
quarkonia arise from the lowest-lying potential; 
hybrid quarkonium states emerge from the ex- 
cited potentials. The leading Schrodinger kinetic 
energy term and the omission of retarded trans- 
verse gluons (in Coulomb gauge), radiated and 
reabsorbed by the slowly moving quarks, could 
be corrected in a systematic expansion. 

The starting point of the second method is 
the lattice formulation of the euclidean NRQCD 
framework with the heavy quark propagating ac- 
cording to a spin-independent nonrelativistic ac- 
tion . Retarded gluon effects are included, spin 
effects and corrections to the leading kinetic term 
can be also added without difficulty. I will report 
new NRQCD results which show that gluon retar- 
dation effects are very small in the bottomonium 
system and its hybrid excitations confirming the 
success of the leading Born-Oppenheimer approx- 
imation. 

The third approach which was explored by 
the MILC collaboration j|] applies Wilson and 
clover lattice fermions in the mass range of the 
charmed quark to explore the hybrid excitations 
of the charmonium system. The heavy c quark 
is treated without nonrelativistic approximation 
and all spin and retardation effects are included. 

I should note that the phenomenology of heavy 
hybrids is very interesting. Early results from the 
CUSB and CLEO collaborations §0] revealed 
a complex resonance structure between the bb 
threshold and 11.2 GeV in e"''e~ annihilation as 
shown in Fig.|l|. This is precisely the energy range 
where the lowest hybrid excitations are expected. 
Today it would only take a very short run for 
the upgraded CLEO detector to rescan this re- 
gion with much higher resolution to find perhaps, 
for the first time, gluon excitations in the hadron 
spectrum under experimentally and theoretically 
rather well-controlled circumstances. This is an 
exciting area where lattice work will remain the 
dominant method. 
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Figure 1. The CUSB scan of the energy region 
above the bb threshold is shown |^ together with 
the expected mass range of the lowest heavy hy- 
brid bbgluon states. The conventional interpre- 
tation of the resonance peaks corresponds to the 
radial excitations of the T particle . 

Born-Oppenheimer Approximation 

The lowest hybrid potential which is de- 
picted in Fig.^ corresponds to a gluon excitation 
with one unit of angular momentum projected 
along the quark-antiquark (molecular) axis with 
CP = — 1 quantum numbers (excited !!„ state 
above the Sg ground state in spectroscopic no- 
tation). I will discuss the full spectrum of gluon 
excitations in Section 3. Following Ref. the 
physical scale rp in Fig.||is defined by the relation 
[r2dV^+(r)/dr] =1.65. 

An important recent study shows in Fig.^ 
that dynamical quark loop effects are not visible 
at larger distances, out to 1.2 fm separation of 
the quark-antiquark pair in the S+ and !!„ po- 
tentials. This further supports the accuracy and 
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Figure 2. Results from mean- field improved 
anisotropic lattice actions ||l^,|l3|] are compared to 
recent (GB) and earlier llsj^PM) simulations 
with Wilson actions. The excited hybrid poten- 
tial has quantum numbers F = IIu. Consistency 
across Refs. 
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Figure 3. Wilson quark loop effects are shown for 
the light quark mass range where m^r ~ mp/2 [ p7[ . 

phenomenological relevance of the approach. 

The spectra of heavy hybrid states built on the 
two lowest hybrid potentials Vn„ (?) and Vj.- (?) 
are shown in Fig.0. The most prominent fea- 



ture of the spectrum is the dense set of radial 
excitations in the rather shallow hybrid poten- 
tials with level separations ~ 200 — 300 MeV. 
The hybrid wavefunctions all vanish at the origin 
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Figure 4. The Schrodinger spectrum of the adi- 
abatic hybrid potentials is shown in comparison 
with the conventional S and P states of the bb 
spectrum. 

due to the repulsive Coulomb core in the poten- 
tial and the special centrifugal term in the adi- 
abatic Schrodinger equation The hybrid 

Schrodinger states are also much more extended 
than the conventional quark-antiquark states, as 
depicted in Fig.^. 

Once the hybrid potentials are determined, 
the Born-Oppenheimer approach yields the en- 
tire leading-order spectrum, in contrast to di- 
rect NRQCD simulations. However, the valid- 
ity of the Born-Oppcnhcimcr approximation re- 
lies on the smallness of retardation effects. In 
addition, spin-dependent corrections have to be 
included for precise predictions. Quantitative 
estimates on retardation effects comes from di- 
rect comparison of mass splittings in the leading 
Born-Oppenheimer approximation with those de- 
termined from NRQCD simulations. 

NRQCD Hybrid Spectrum 

Recent results on heavy hybrid states were 
reported using anisotropic coarse lattices with 
mean- field improved action |l3,n3|, and Wilson 
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Figure 5. The radial probability distribution in 
the lowest IIu hybrid state is shown in compar- 
ison with conventional bb S and P states of the 
bottomonium (upper box) and charmonium sys- 
tem. 

gauge action |]l8|,^. In Ref. the bare quark 
mass was taken to be agMb = 2.56. The so-called 
kinetic mass of the T was then determined from 
its low-momentum dispersion relation. Half of 
this mass was used for the quark mass in the 
leading Born-Oppenheimer calculation. This en- 
sured that the T kinetic masses were identical in 
both calculations. A typical effective mass plot 
for hybrid states from the NRCQD run with im- 
proved action on asymmetric lattices is shown in 
Fig.^ fl^. The hybrid mass with Wilson gauge 
action |18| was found noticeably higher, but part 
of the disagreement might be attributable to the 
different choice of setting the physical scale. Re- 
sults from Ref. |l3j are consistent with earlier 
work [ p"2|jTg| ] adding improved control on system- 
atics and statistics. One of the new features here 
is the clear signal for the first radial hybrid ex- 
citation. The statistical error on the mass of the 
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Figure 6. Mass estimates are shown for the lowest 
hybrid IIu state and its first radial excitation jl^ . 



hybrid bbg state is small but the systematic study 
of lattice artifacts remains incomplete. Assuming 
small lattice artifacts in the simulation results, 
it was shown that gluon retardation affects the 
spin-averaged mass splittings by less than ten per- 
cent, validating the leading Born-Oppenheimer 
approximation |l3| . 

Results from the MILC Collaboration 

The MILC collaboration repeated its original 
Wilson simulations of the hybrid meson spec- 
trum 1^ using the clover action The new 
results, shown in Fig.^ for the 1 heavy hybrid 
mass, are consistent with earlier MILC Wilson 
runs providing an important check on charmed 
quark lattice artifacts. Preliminary results for the 
wave function of the 1 state in Coulomb gauge 
were also reported at the conference S. 



2. GLUEBALL SPECTRUM 

In the perspective of ten years it is quite re- 
markable how much is known now about the 
quenched glueball spectrum. Hopefully, it will 
not take another ten years before quark-antiquark 
vacuum polarization effects are included system- 
atically, with control comparable to the quenched 
results. 
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Figure 7. The lowest hybrid states are shown 
in the charmonium and bottomonium spec- 
trum from three different methods described in 
the text. The spin-dependent potential model 
spectrum is superimposed for comparison from 
Rcf. where the notation is explained. 

Review of Recent Lattice Results 

The construction of improved anisotropic lat- 
tice actions using mean-field improved perturba- 
tion theory is discussed in Ref. with further 
pointers to the literature. All recent results re- 
viewed here are based on this action which was 
designed to be 0{al) improved and to include 
terms that directly couple links only on adja- 
cent time slices to avoid spurious modes in the 
gluon propagator. The hadronic scale param- 
eter ro, defined earlier, is used in the conver- 
sion of the measured glueball masses into phys- 
ical units. The method has been applied to 
the SU(3) glueball spectrum by Morningstar and 
Peardon p2| w hose comprehensive lattice spin 
analysis ||23| , ^4| is shown in Fig|8[ Recently a de- 
tailed analysis of the SU(2) glueball spectrum us- 
ing very similar methods was also presented psf . 
The assignment of continuum spin labels to lat- 
tice states in Fig.|| is based on transformation 
properties of the glueball states under the cubic 
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Figure 8. The SU(3) lattice spectrum of glueball 
states in the quenched approximation. 



point group, combined with parity and charge 
conjugation operations [ p2|j2^ ]. The continuum 
spin quantum numbers are not always fully re- 
solved across cubic group representations (des- 
ignated as ambiguous states). The bottom of 
the shaded area in each channel of Fig.|| indi- 
cates where the continuum threshold of the light- 
est two-glueball state is located. Glueball states 
above the continuum threshold are not included 
in the final spin analysis of the continuum spec- 
trum in Figj^. The two lightest glueballs, the 
scalar and tensor, have candidate resonances p7| - 
!§: the /o(1500) and /j(1710) scalar states and 
the ^(2230) tensor candidate. Both /o(1500) and 
/j(1710) are consistent with lattice mass predic- 
tions in the pure-gauge theory |^,|^|2| of 1600 
McV with systematic errors of approximately 100 
MeV. The glueball spectrum from the Wilson ac- 
tion is shown in Fig.|l0|, together with the spec- 
trum from a constituent gluon model. 

The SESAM collaboration Q has investigated 
the scalar and tensor glueballs on their large en- 
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Figure 9. Final spin analysis of the continuum 
glueball spectrum. 



12 



10 



SU(3) Glueball Spectrum 

Wilson action data and gluon modes 



(TM)3 I 



|(TM)(TMi 2 I 



tTMXTM) |(TE)(TM)| f 



- 1^- ■! 

■ g- i (™)3 r 

(TE* I 



l(TE)(TE)| 



J^'^ continuum limit extrapolations 



4000 



3000 



2000 



1000 



+ +- 

PC 



Figure 10. The range of energy levels from several 
Wilson action runs is indicated by shaded areas. 
The constituent gluon spectrum (discussed in the 
text) is also displayed with the occupied gluon 
modes. 



semble of unquenched (Wilson gauge and quark 
action) configurations. Only little change has 
been found in the masses of these states as the 
sea quark mass runs down from the strange quark 
mass to (m^/TOp w 0.5) but enhanced finite- 
volume dependence was noticed which is perhaps 
attributable to dynamical quark effects. 

Composite Glueball Operators 

In a confining theory it is difficult to relate 
physical glueball states to fundamental gluon 
fields. There is, however, a qualitative expecta- 
tion that higher dimensional interpolating opera- 
tors will create higher mass states p^-|3^]. This 
rule works well for light mesons and baryons and 
now I will apply it to the glueball spectrum. 

Local composite operators of the gluon field act 
as interpolating fields for glueball states. For ex- 
ample, the state with quantum numbers is 
created from the physical vacuum by the opera- 
tor Ea ■ Ba (a=l,2,...,8 stands for color summa- 
tion) so that (0 ^|Ea ■ Ba|vacuum) ^ 0. It is 
not difficult to determine the quantum numbers 
of all the glueball states which are generated by 
the lowest dimensional (dim=4,5,6) operators. A 
small representative sample is listed in Table |^. 

Table 1 

Sample of dim=4,6 operators which create vari- 
ous J^*-" quantum numbers. 



dimension 




Operator 


4 


0-+ 


Ea • Ba 


4 


0++ 


E^-B2 


4 


2++ 




6 


0-+ 


fabc(Ea X Eb) • Ec 


6 


1-+ 


fabc(Ba • Eb)Bc 


6 


2-+ 


fabc(Ea X Eb)^E^e + - 



In classifying the lowest glueball excitations 
according their quantum numbers, Morningstar 
and Peardon report 15-1-5 states where the last 
five are labelled as questionable. The lowest di- 
mensional operators with dim=4,5,6 create states 
whose quantum numbers match onto the observed 
spectrum nearly perfectly. This is quite remark- 
able, given the simplicity of the argument. 



Constituent Gluon Model 

The notion of consitucnt gluons stems from 
the general idea of quasi-particles in some ef- 
fective QCD mean field theory of hadrons. In 
such a description the fundamental gluon field is 
replaced by Hartree modes of some constituent 
gluon field with residual perturbative interactions 
inside glueballs. It is hoped, for example, that the 
1 /N expansion might eventually lead to a similar 
Hartree picture. I will identify the Hartree states 
with free gluon modes in a bubble (bag) with con- 
fining boundary conditions ||3^,^. Admittedly 
speculative, the model in its simplest form leads 
to detailed predictions on the glueball spectrum. 

We will picture the QCD vacuum as a diaelec- 
tric medium (dual superconductor) with vanish- 
ing chromoelectric permeability and infinite mag- 
netic permeability. Glueballs arc fioating bubbles 
in the physical vacuum (Fig.p^) and their per- 
turbative interior is occupied by constituent TE 
and TM gluon modes appropriate for the con- 
fining boundary conditions of a dual supercon- 
ductor mjl^. The energy density in the physi- 
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Figure 11. Hartree states of TE and TM gluon 
modes. 



cal vacuum (which is the required vacuum energy 
to create the bubble) is given by the well-known 
expression B = — 3|;fas(0|F^j,|0) from the trace 
anomaly of the energy-momentum tensor [^7| . 
The quantum numbers and energies of the lowest 
gluon modes are given in Table || and the fields 
of the TE and TM modes are depicted in Fig, 11. 
The energy of the localized glueball is calculated 

Table 2 

Lowest TE and TM gluon modes in spherical bub- 
ble with energies Xi/R for the i'^ mode. 



1 


mode 




Xi 


1 


TE 


1+ 


2.744 


2 


TE 


2- 


3.96 


1 


TM 




4.49 



from 

Egiucbaii = yR=^B+^niXi/R-c^AtA^SiS2, (1) 

i 

where Uj is the occupation number of the i*'' gluon 
mode and the numerical constant c is close to 
1/4. The last term in Eq. (|l|) represents the 
interaction energy of the gluon modes. I chose 
B^/** — 280 MeV for the vacuum energy den- 
sity which is appropriate for heavy quark spec- 
troscopy and the hybrid potentials Q (B^/^ ~ 
250 MeV — 350 MeV is the preferred range from 
the evaluation of the gluon condensate |3^ ) . The 
mass of the glueball has to be corrected for the 
localization of the center-of-mass, 

Mglueball = Egi^jji^jjjj — ^ Uixf /R^, (2) 
i 

where the second term is an estimate of the ki- 
netic energy (P^). The energies and radii of the 
various glueball states are determined from the 
minimization of Eqs.(0j|) with respect to R. The 
spectrum is displayed in Fig.|l^ with as(l/R) = 
0.5 where R ~ 0.5 fm. 

The old bag calculations of the glueball spec- 
trum did not work. They led to very low glue- 
ball masses with wrong spin splittings for two 
reasons: the wrong value of B and an unaccept- 
ably large effective as were used simultaneously 
from light hadron spectroscopy Jss],^. I have 



nothing to say about light hadrons which are in- 
timately connected with the chiral condensate of 
the vacuum. In contrast, the dual superconduc- 
tor picture gives consistent phenomenology for so 
many aspects of the microscopic vacuum struc- 
ture in the gluon sector of the theory |^^. It 
remains a mystery why this simple caricature of 
glueball states works so well in describing the 
quenched spectrum. Only a microscopic theory 
of the vacumm could shed light on this question. 

3. THE CONFINING FLUX 

I will discuss first an interesting microscopic 
vacuum model to illustrate the generic features 
of the confining flux. The Lagrangian of the 
same model was also used in two contributions to 
the conference in the study of flux flssion (string 
breaking) plEol. 



Microscopic Vacuum Model 

We will investigate the response of the micro- 
scopic monopole vacuum to Wilson loops in the 
2-1-1 dimensional Georgi-Glashow model as the 
simplest example of quark confinement pl[ |. In 
the Higg vacuum of the theory the SU(2) gauge 
symmetry is broken to U(l) with a neutral and 
massless gauge boson A^. The massive Higgs bo- 
son is frozen out from the low energy theory, to- 
gether, in the weak coupling limit g^/mw ^ 1, 
with the two massive gauge bosons W^ where 
the dimensional gauge coupling g sets the phys- 
ical scale. The euclidean vacuum is dominated 
by a dilute gas of instantons which are 't Hooft- 
Polyakov monopoles interacting with a magnetic 
Coulomb force. 

The monopole vacuum responds to external 
U(l) charges by flux formation. This is calculable 
from the response of the vacuum to an external 
magnetic potential 77, given by the partition func- 
tion 



Zirj) = / Pxcxp - 



167r2 
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d'^r 



(3) 



in the sine-Gordon field representation of the 
monopole plasma |41[. The sine-Gordon mass 



M is expressed in terms of the fugacity C of 
the monopole plasma by — 327r2(^/g2 where 
^ ~ exp(-— const • mw/g^) is exponentially small 
in the weak coupling limit so that we have 
a dilute monopole system. A pair of static 
sources corresponds to an external magnetic po- 
tential rywi which is generated by an imagi- 
nary current I = i • g/2 running around the Wil- 
son loop immersed in the monopole plasma, 
(exp(i / A^dx^)) = Z(?7wi)/Z(0) . In the weak 
coupling limit the semiclassical expansion is ap- 
plicable (Debye-Hiickel theory), and flux forma- 
tion in leading order is described by the solution 
of the sine-Gordon equation in the presence of 
the current loop. An exact sine-Gordon soliton 
solution at infinite separation of the sources, 



Xsoliton 
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(4) 



describes a confining fiux running along the x 
direction, with the transverse profile given by 
Xsoliton in the y-direction. This solution in 2-1-1 di- 
mensional Minkowski space-time corresponds to 
the electric flux profile 
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where the width of the flux which contains 
more than 90 percent of the flux energy is ap- 
proximately 4M, set by the Debye-Hiickel (sine- 
Gordon) mass M. The asymptotic value of the 

2 

string tension is given by cr = ^f^^- 

At finite separation of the sources we solve the 
sine-Gordon equation numerically with typi- 
cal solutions depicted in Fig.|l^. We find a grow- 
ing and stretching bubble in the vacuum crossing 
from oblate to prolate shape and asymptotically 
converging to the shape and intrinsic profile of the 
exact solution of Eq. (||) between the two fixed 
sources but not close to them. The vortex energy 
is shown in Fig.|l^. When the separation of the 
sources becomes comparable to the width of the 
flux, the rapid onset of linear behavior in the flux 
energy, well before the asymptotic form of the flux 
is reached, is a remarkable feature of Figs. 12 1^. 
It is likely that a rather similar situation occurs in 
QCD where we would like to answer the question: 
What sets the scale of string formation? 



Flux formation in monopole vacuum 
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Figure 12. Flux profiles at small and intermediate 
separation of the sources. The bottom part is a 
contour plot of distortions on the gutter shaped 
effective Schrodinger potential as discussed in the 
text. 




Figure 13. Flux energy with arrows indicating 
the flux configurations of Fig.|2|. The slope of 
the dashed line is the exact string tension of the 
infinitely long flux. 



into the details of the excitation spectrum. One 
can excite the color sources at the two ends by 
putting the quark and the antiquark into higher 
representations of the SU(3) color group, or inves- 
tigate flux formation with several quark sources. 

Casimir scaling was a much discussed topic at 
the conference: according to folklore the string 
tension should be proportional to the value of 
the Casimir operator in the color group repre- 
sentation. Deldar reported results for the flux 
energy in the sextet and octet color representa- 
tions of the static quark-antiquark pair |^^. A 
linearly rising potential with Casimir scaling was 
observed in the approximate range of 0.5 fm to 
1.2 fm source separation, even for octet sources 
where gluon screening is expected to split the flux 
asymptotically. This is in agreement with earlier 
findings [Hg. 

Although it is popular to interpret Casimir 
scaling at intermediate QQ-separation as a test of 
microscopic confinement mechanisms (Z(N) flux, 
monopoles, etc.), the test remains incomplete 
without a more detailed microscopic understand- 
ing of flux formation and its dependence on color 
representations |46|. 



In some interesting work in progress |21|, the 
potential energy of the 3Q system was studied. 
Anticipating the formation of a Y-shape string at 
large separation, the potential is plotted in Fig.|lJ 
as a function of inter-quark separation when the 
three quark sources are located in a plane at equal 
angles. The lesson from the growing and stretch- 
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Casimir Scaling and Multiple Sources 

Several interesting tests can be done in QCD 
on the flux formation mechanism without looking 



Figure 14. 3Q potential in Y-configuration. 



ing bubble is that the rapid onset of the hnear 
potential cannot be interpreted as evidence for 
the early formation of a Y-shape string. Earlier 
work in the bag model illustrates that a very sim- 
ilar potential energy shape develops precociously 
before string formation H] , as depicted in Fig.[l5|. 
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Figure 15. 3Q bag potential in Y-configuration. 
The dashed line indicates spherical approxima- 
tion. 

Interesting new results were also reported at 
the conference on 4Q flux configurations 
which hopefully will advance our understanding 
of flux formation in the confining vacuum. 

Where Is the String Limit? 

The confining flux is represented by the quan- 
tum sine-Gordon soliton in the monopole vac- 
uum. Excitations of the flux are given by the 
spectrum of the fluctuation operator M = — V^-f 
U" (Xsoiiton) where U(x) is the field potential en- 
ergy of the sine-Gordon field. The spectrum of 
the fluctuation operator Ai is determined from 
a two-dimensional Schrodinger equation at fi- 
nite separation of the sources. In the limit of 
asymptotically large rgq, the equation becomes 
separable in the longitudinal x and transverse 
y coordinates. The transverse spectrum is in 
close analogy with the quantization of the one- 
dimensional classical sine-Gordon soliton. There 
is always a discrete zero mode in the spectrum 



which is enforced by translational invariance in 
the transverse direction. Scattering states of mas- 
sive gauge bosons ( "glueballs" ) on the confining 
flux are separated by the mass gap M of the sine- 
Gordon theory. 

Plane wave excitations with momentum k along 
the longitudinal direction, when combined with 
the zero mode of the transverse spectrum, cor- 
respond to Goldstone excitations along the in- 
finitely long flux. This spectrum is mapped into 
the spectrum of a string in the semiclassical quan- 
tization procedure |49|. In string theory, the 



string with fixed ends is described by Dirichlet 
boundary conditions which would correspond to 
closing down on the two ends the long effective 
Schrodinger potential well with two infinite po- 
tential walls (a gutter between infinite walls). 

At finite Tqq, this picture, however is grossly 
distorted. Far from the sources and from the 
main axis, the effective Schrodinger potential well 
always has a finite height M. Therefore the gutter 
should be closed by two walls of finite height M 
which are penetrable and distort the Goldstone 
spectrum. In addition, as depicted in Fig.|l^, the 
potential well is further deformed close to the 
sources. As a consequence, at smaller separation 
of the sources, the fluctuation spectrum is not 
string-like at all. At large separations, however, 
the fluctuation spectrum of the quantum soliton 
(confining flux) will smoothly deform with grow- 
ing rQQ into the expected Goldstone spectrum. 
Based on the above picture this is not expected 
before the length of the flux exceeds several times 
the width of the flux. 

In QCD the robust features of the Goldstone 
spectrum of the confining flux are described by an 
effective string theory where the two-dimensional 
vector ^ measures the deviations of the flux from 
its straight line position |50|. The effective low 



energy expansion is based on the assumption that 
massive intrinsic excitations of the flux are decou- 
pled at low energies from the low-energy string- 
like excitations. The effective Lagrangian of this 
description has two important symmetries: Eu- 
clidean invariance in the world sheet variables 
and 0(2) translations and rotations in the ^ field 
variable. The second symmetry requires that 
CcsiO only depend on the derivatives of ^, 



CeS = ^ctdfj,£^ ■ + where the dots represent 
terms with four, or more derivatives and a is a 
dimensional constant. For low frequency excita- 
tions the first term in Ccs{S,) dominates all the 
higher dimensional operators. 

The Gaussian effective Lagrangian of the fiuc- 
tuating fiux describes massless Goldstone excita- 
tions which are associated with the restoration of 
Euclidean translations and rotations due to the 
formation of the confining quantum flux. The 
full non-linear realization of euclidean symmetries 
on ^ would amount to the summation of higher 
derivative operators requiring special coefficients 
and resulting in the Nambu-Goto action, 

-Coff = a\J 1 + df,^- df,^= ^a9^f • d^^+ ... (5) 

This is in analogy with the nonlinear chiral La- 
grangian of massless pions. However, the spe- 
cial higher derivative terms in the expansion 
of the Nambu-Goto action will compete with 
other higher dimensional operators (rigidity, tor- 
sion, and some other properties of the confining 
fiux). It is, therefore, difficult to believe that the 
Nambu-Goto action would play any useful role in 
the description of the excitation spectrum of the 
confining flux. 

The ground state energy of a string of length L 
with clamped ends is given by 

Eo(L)=f7.L-^ + ... (6) 

in 1-loop approximation where a designates the 
string tension and the second term represents 
the Casimir energy of zero-point fluctuations [^0| . 
The excitation energies above the ground state 
have the simple string spectrum, 

E„ = Eo + n-^, n=l,2,3,... (7) 

It was repeatedly claimed by lattice practition- 
ers that the Casimir term was seen in the static 
quark-antiquark potential (string ground state) 
for less than two fermi separation, interpreted as 
strong evidence for string formation in QCD. I 
will show now that this claim was intrinsically 
flawed. 

In QCD we do not have the underlying mi- 
croscopic theoretical picture of the confining vac- 



uum, but we expect to learn about it by the deter- 
mination of the excitation spectrum of the con- 
fining fiux in lattice simulations. The first step 
in determining the spectrum of the stationary 
states of the gluon field in the presence of the 
static quark-antiquark pair, fixed in space some 
distance apart, is to classify the energy levels 
in terms of the symmetries of the problem [p^ . 
Such a system has cylindrical symmetry about 
the axis f passing through the static quark and 
the antiquark (the molecular axis). The total an- 
gular momentum of the gluons is not a conserved 
quantity, but its projection onto the molecular 
axis is and can be used to label the energy levels 
of the gluon field. We adopt the standard nota- 
tion from the physics of diatomic molecules and 
use A to denote the magnitude of the eigenvalue 
of the projection Jg- r of the total angular mo- 
mentum Jg of the gluon field onto the molecular 
axis f . The capital Greek letters E, 11, A, $, . . . 
are used to indicate states with A = 0, 1, 2, 3, . . ., 
respectively. The energy of the gluon field is un- 
affected by refiections in a plane containing the 
molecular axis; since such refiections interchange 
states of opposite handedness, given by the sign 
of Jg - i, such states must necessarily be degen- 
erate (A doubling). However, this doubling does 
not apply to the S states; E states which are even 
(odd) under a refiection in a plane containing the 
molecular axis are denoted by a superscript -t- 
(— ). The combined operations of charge conju- 
gation and spatial inversion about the midpoint 
between the quark and the antiquark is also a 
symmetry and its eigenvalue is denoted by rjcp. 
States with ry^p = 1(— 1) are denoted by the sub- 
scripts g (u). Hence, the low- lying levels are la- 
belled S+, , S+, Ug, H„, Ag, A„, and 
so on. For convenience, we use F to denote these 
labels in general. 

The continuum-limit extrapolations are 
shown in Fig.|l^. The ground state is the 
familiar static-quark potential. A linearly-rising 
behaviour dominates the potential once r ex- 
ceeds about 0.5 fm and we find no deviations from 
the linear form up to 4 fm. The lowest-lying ex- 
citation is the Hu. There is definite evidence of a 
band structure at large separation: the S+', Ug, 




Figure 16. The excitation spectrum of the con- 
fining flux. The energy levels of N=l,2,3,4 Gold- 
stone excitations are shown in the insert. 



and Ag form the first band above the n„; the E+, 
E~, n^/<i>„, and A„ form another band. The S~ 
is the highest level at large r^^. This band struc- 
ture breaks down as the separation of the sources 
decreases below 2 fm. In particular, two levels, 
the and E", drop far below their degenerate 
partners as the separation between sources be- 
comes small. Note that for above 0.5 fm, all 
of the excitations shown are stable with respect to 
glueball decay. As the separation of the sources 
decreases below 0.5 fm, the excited levels even- 
tually become unstable as their gaps above the 
ground state S+ exceed the mass of the lightest 
glueball. 

The level orderings and approximate degen- 
eracies of the gluon energies at large separation 
match those expected of the Goldstone modes. 
However, the precise Goldstone gap behaviour is 
not observed. The two S~ states are in violent 
disagreement with expectations from a fluctuat- 
ing string. Note also that the results clearly dis- 
agree with the energy spectrum of a Nambu-Goto 
string naively applied in four continuous space- 
time dimensions. 

These results are rather surprising and cast se- 
rious doubts on the validity of treating glue in 
terms of a fluctuating string for quark-antiquark 



separations less than 2 fm. Note that such a con- 
clusion does not contradict the fact that the E+ 
ground state energy rises linearly for rgq as small 
as 0.5 fm. A linearly-rising term is not necessar- 
ily indicative of a string as we have seen in earlier 
examples. For t^qq greater than 2 fm, there are 
some tantalizing signatures of Goldstone mode 
formation, yet significant disagreements still re- 
main. To what degree these discrepancies can be 
explained in terms of distortions of the Goldstone 
mode spectrum arising from the spatial fixation 
of the quark and antiquark sources (clamping ef- 
fect) is currently under investigation ji^ ]. 

For reasons explained here I remain puzzled 
and reserved on reported results that large Wil- 
son loops in lower dimensional models agree with 
the predictions of a fluctuating string |5^ even 
for relatively small separation of the sources. 

Flux Fission (String Breaking) 

The confining flux is expected to fission into a 
pair of static mesons, also known as string break- 
ing, when its energy is large enough to pair pro- 
duce light dynamical quarks. At large separation 
of the color sources the static potential, as cal- 
culated from Wilson loops, describes the force 
between static mesons which are color sources 
screened by light quark fields (like the static B 
meson in the infinite b-quark mass limit). We 
expect the linearly rising potential to asymptot- 
ically cross over into Yukawa form controlled by 
the lowest mass exchange (pion). The attractive, 
or repulsive nature of the Yukawa force will de- 
pend on the spin-isospin quantum numbers of the 
static BB pair which is not taken into account in 
Wilson loop operators. The crossover range is 
expected to occur at a separation where the en- 
ergy of the confining flux exceeds twice the static 
meson energy. 

In QCD simulations with two flavors of dy- 
namical quarks, where the crossover range has 
been reached, new results were reported at the 
conference (Fig jl7|) w ithout any visible string 
breaking effects 51 5|] . It was suggested us- 
ing strong coupling ideas [|3| that string break- 
ing is a mixing phenomenon, involving both the 
string and the static two-meson state. Thus, in 
order to confirm the mixing picture, the conven- 
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Figure 17. Dynamical fcmiion loop effects on the 
static potential. 

tionally used Wilson loops have to be supple- 
mented by explicit two-meson operators. This 
idea was illustrated in two contributions to the 
conference studying the Georgi-Glashow 

model with the Higgs field in the fundamental 
representation playing the role of the screening 
matter field. In Fig.|l^ results from the three- 
dimensional model are shown [^9| with findings 
very similar to the four-dimensional study [^0| . 
Since the quenched static potential between ad- 
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Figure 18. Mixing effects between the Wilson 
loop operator and two static mesons. 

joint sources is screened by the gluon field, its 
study would also provide further tests on string 



breaking 

It was reported at the conference from the 
study of Polyakov loop correlators, where mixing 
with the static meson pair is different from Wilson 
loops, that the finite temperature static potential 
does exhibit screening behavior |5^. Since the 
calculation was done just below the deconfining 
temperature it would be interesting to extend this 
work to lower temperatures. It was also reported 
that the screening of large Wilson loops was easy 
to detect in lower dimensional models |55 5^. 
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